Observation of a new charmonium-like state produced in association 
with a J/ip in e + e annihilation at ^/s « 10.6 GeV 

K. Abe, 7 I. Adachi, 7 H. Aihara, 40 K. Arinstein, 1 Y. Asano, 44 V. Aulchenko, 1 T. Aushev, 11 T. Aziz, 36 
A. M. Bakich, 35 V. Balagura, 11 M. Barbero, 6 I. Bedny, 1 U. Bitenc, 12 I. Bizjak, 12 A. Bondar, 1 M. Bracko, 7 ' 18 > 12 
J. Brodzicka, 25 T. E. Browder, 6 Y. Chao, 24 A. Chen, 22 B. G. Cheon, 2 R. Chistov, 11 S.-K. Choi, 5 Y. Choi, 34 
Y. K. Choi, 34 A. Chuvikov, 30 S. Cole, 35 J. Dalseno, 19 M. Danilov, 11 M. Dash, 45 A. Drutskoy, 3 S. Eidelman, 1 
D. Epifanov, 1 S. Fratina, 12 N. Gabyshev, 1 T. Gershon, 7 G. Gokhroo, 36 B. Golob, 17 - 12 H. C. Ha, 14 J. Haba, 7 
Y. Hasegawa, 33 K. Hayasaka, 20 H. Hayashii, 21 M. Hazumi, 7 L. Hinz, 16 Y. Hoshi, 38 S. Hou, 22 W.-S. Hou, 24 
Y. B. Hsiung, 24 T. Iijima, 20 A. Ishikawa, 7 M. Iwasaki, 40 Y. Iwasaki, 7 P. Kapusta, 25 T. Kawasaki, 27 H. R. Khan, 41 
H. J. Kim, 15 S. M. Kim, 34 K. Kinoshita, 3 S. Korpar, 18 - 12 P. Krokovny, 1 C. C. Kuo, 22 A. Kuzmin, 1 J. S. Lange, 4 
G. Leder, 10 T. Lesiak, 25 A. Limosani, 7 S.-W. Lin, 24 D. Liventsev, 11 F. Mandl, 10 T. Matsumoto, 42 

A. Matyja, 25 Y. Mikami, 39 W. Mitaroff, 10 H. Miyata, 27 R. Mizuk, 11 Y. Nagasaka, 8 E. Nakano, 28 M. Nakao, 7 
S. Nishida, 7 O. Nitoh, 43 S. Ogawa, 37 T. Ohshima, 20 S. Okuno, 13 S. L. Olsen, 6 H. Ozaki, 7 P. Pakhlov, 11 

H. Palka, 25 C. W. Park, 34 H. Park, 15 K. S. Park, 34 L. S. Peak, 35 L. E. Piilonen, 45 A. Poluektov, 1 Y. Sakai, 7 
N. Sato, 20 T. Schietinger, 16 O. Schneider, 16 A. J. Schwartz, 3 R. Seidl, 31 K. Senyo, 20 M. E. Sevior, 19 

B. Shwartz, 1 V. Sidorov, 1 A. Somov, 3 R. Stamen, 7 M. Staric, 12 T. Sumiyoshi, 42 S. Y. Suzuki, 7 F. Takasaki, 7 
N. Tamura, 27 M. Tanaka, 7 Y. Teramoto, 28 X. C. Tian, 29 K. Trabelsi, 6 T. Tsukamoto, 7 S. Uehara, 7 T. Uglov, 11 

K. Ueno, 24 Y. Unno, 7 S. Uno, 7 P. Urquijo, 19 G. Varner, 6 K. E. Varvell, 35 S. Villa, 16 C. C. Wang, 24 
C. H. Wang, 23 M.-Z. Wang, 24 E. Won, 14 Q. L. Xie, 9 B. D. Yabsley, 45 A. Yamaguchi, 39 Y. Yamashita, 26 
M. Yamauchi, 7 J. Ying, 29 C. C. Zhang, 9 J. Zhang, 7 L. M. Zhang, 32 Z. P. Zhang, 32 and V. Zhilich 1 

(The Belle Collaboration) 

1 Budker Institute of Nuclear Physics, Novosibirsk 

2 Chonnam National University, Kwangju 
University of Cincinnati, Cincinnati, Ohio 45221 
' University of Frankfurt, Frankfurt 
5 Gyeongsang National University, Chinju 
6 University of Hawaii, Honolulu, Hawaii 96822 
7 High Energy Accelerator Research Organization (KEK), Tsukuba 
8 Hiroshima Institute of Technology, Hiroshima 
9 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 
10 Institute of High Energy Physics, Vienna 
11 Institute for Theoretical and Experimental Physics, Moscow 
12 J. Stefan Institute, Ljubljana 
13 Kanagawa University, Yokohama 
14 Korea University, Seoul 
15 Kyungpook National University, Taegu 
16 Swiss Federal Institute of Technology of Lausanne, EPFL, Lausanne 
17 University of Ljubljana, Ljubljana 

18 University of Maribor, Maribor 
19 University of Melbourne, Victoria 

20 Nagoya University, Nagoya 
21 Nara Women's University, Nara 
22 National Central University, Chung-li 
23 National United University, Miao Li 
^Department of Physics, National Taiwan University, Taipei 
H. Niewodniczanski Institute of Nuclear Physics, Krakow 
2b Nippon Dental University, Niigata 
21 Niigata University, Niigata 
28 Osaka City University, Osaka 
29 Peking University, Beijing 
80 Princeton University, Princeton, New Jersey 08544 
'"RIKEN BNL Research Center, Upton, New York 11973 
32 University of Science and Technology of China, Hefei 
Shinshu University, Nagano 



2 



3J> Sungkyunkwan University, Suwon 
35 University of Sydney, Sydney NSW 
36 Tata Institute of Fundamental Research, Bombay 
s7 Toho University, Funabashi 
38 Tohoku Gakuin University, Tagajo 
39 Tohoku University, Sendai 
■*° Department of Physics, University of Tokyo, Tokyo 
J > 1 Tokyo Institute of Technology, Tokyo 
,fi 'Tokyo Metropolitan University, Tokyo 
* Tokyo University of Agriculture and Technology, Tokyo 
^ University of Tsukuba, Tsukuba 
* Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 

We report the first observation of a new charmonium-iike state at a mass of (3.943 ± 0.006 ± 
0.006) GeV/c 2 . This state, which we denote as X(3940), is observed in the spectrum of masses 
recoiling from the J/ip in the inclusive process e + e~ — > J/ip + anything. We also observe its decay 
into D* D and determine its intrinsic width to be less 52MeV/c 2 at the 90% C.L. These results are 
obtained from a 357 fb" 1 data sample collected with the Belle detector near the T(AS) resonance, 
at the KEKB asymmetric-energy e + e~ collider. 

PACS numbers: 13.66.Bc,12.38.Bx,14.40.Gx 



Recently there have been a number of reports of 
new charmonium or charmonium-like states: rj c (2S) 0, 
A(3872) 0, y(3940) and F(4260) 0. The latter 
three states have not been assigned to any charmonium 
states in the conventional quark model. Moreover, char- 
monium production in different processes is not well 
understood. One striking example is the surprisingly 
large cross section for double charmonium production in 
e + e~ annihilation observed by Belle 5] and confirmed 
by BaBar . These experimental results have generated 
renewed theoretical interest in the spectroscopy, decays 
and production of charmonium. 

In this paper we report the observation of a new 
charmonium-like state above DD threshold, A (3940), 
produced in the process e + e~ — > J / ijj X (39 AO) . We also 
present results from searches for A" (3940) decay into DD, 
D*D and J/ipoj. The data used for this analysis corre- 
spond to an integrated luminosity of 357 fb -1 collected 
by the Belle detector at the T(4S) resonance and nearby 
continuum at the KEKB asymmetric-energy e + e~ col- 
lider. 

The J/ip reconstruction procedure is identical to our 
previously published analyses 0, 0] . Oppositely charged 
tracks that are positively identified as muons or electrons 
are used for J/ip — » £ + £~ reconstruction. A partial cor- 
rection for final state radiation and bremsstrahlung en- 
ergy loss is performed by including the four-momentum 
of every photon detected within a 50 mrad cone around 
the electron direction in the e + e~ invariant mass calcu- 
lation. The two lepton candidate tracks are required to 
have a common vertex, with a distance from the inter- 
action point in the plane perpendicular to the beam axis 
smaller than 1 mm. The J/ip — > l + £~ signal region is de- 
fined by the mass window \M e + e - — Mj/^\ < 30MeV/c 2 
(« 2.5 a). J/ip candidates in the signal window are sub- 
jected to a mass and vertex constrained fit to improve 



their momentum resolution. QED processes are substan- 
tially suppressed by requiring the total charged multi- 
plicity (iV c h) in the event to be greater than four. Back- 
ground due to J/ip mesons produced from BB events is 
removed by requiring the center-of-mass (CM) momen- 
tum P*j/^ to be greater than 2.0GeV/c. As in the previ- 
ous analysis, we define the recoil mass as 

M Iecoil (J/iP) = ^(Ecus-E^-p*^, (1) 

where Ej,^ is the J/ip CM energy after the mass con- 
straint. 

For the study of the A(3940) -> D<~*W, we recon- 
struct D° candidates using five decay modes: K~ir + , 
K~K+, K-TT-TT+TT+, K q s tt+tt- and K-jt+tt - and D+ 
candidates using K~ir + n + , K~K + n + and Kgir + . For 
the D° — > K~7r~ir + 7r + and D° — > K~tt + it modes, mass 
windows of ±10MeV/c 2 and ±20MeV/c 2 are used; a 
±15 MeV/c 2 mass window is used for all other modes (ap- 
proximately 2.5 a in each case) . To improve their momen- 
tum resolution, D candidates are refitted to the nominal 
D° or D + masses. To study the contribution of combina- 
torial background under the D peak, we use D sidebands 
with mass windows that are four times as large. For the 
A (3940) — > J/ipoj search, candidate u> mesons are recon- 
structed from 7r + 7T~7r° combinations within ±20 MeV/c 2 
(~ 2.5 er) of the nominal u> mass. The uo sideband region 
is defined by 30 < \M(tt+it- tt q ) - M u \ < 50MeV/c 2 . 

The recoil mass spectrum for the inclusive J/ip event 
sample is shown in Fig. Here, in addition to the 
three previously reported peaks at the r/ c , Xco an d rj c (2S) 
masses, there is a fourth peak above DD threshold. We 
perform a fit to this spectrum that includes the three 
previously seen charmonium states plus a fourth state 
with mass near 3.94GeV/c 2 . The expected signal line- 
shapes are determined from Monte Carlo (MC) simula- 
tion as described in previous Belle publications 0, 0]. 
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FIG. 1: The distribution of masses recoiling against the recon- 
structed J/ip in inclusive e + e~ — > J/ipX events. The curves 
are described in the text. 

The mass values for the r] c , x c0 , r] c (2S) and X(3940) are 
free parameters in the fit, the widths of r\ c and Xco are 
fixed to PDG values || and the i] c (2S) width is fixed to 
T = 17MeV/c 2 0. The width of the new state is a free 
parameter in the fit. The signal function for the A(3940) 
is a convolution of the (zero-width) MC line shape with 
a Breit-Wigner function. The background is parameter- 
ized by a second order polynomial and a threshold term 
( \ZM icco u( J/ip) — 2Mo) to account for a possible contri- 
bution from e+e- -> J/tpD^D^. 

The fit results are given in Tabled The significance for 
each signal is defined as y— 2 ln(£o/£max), where Co and 
£max denote the likelihoods returned by the fits with the 
signal yield fixed at zero and at the fitted value, respec- 
tively. The significance of the A'(3940) signal is 5.0 a. 
The fitted width of the A (3940) state is consistent with 
zero within its large statistical error: T = 39±26 MeV/c 2 . 
The fit results are shown in Fig. ^ as the solid curve; the 
dashed curve is the background function. 



TABLE I: Summary of the signal yields, charmonium masses 
and significances for e + e~ — s- J/ip (cc) res . 



(cc) rcs 


N 


M [GeV/c 2 ] 


iV CT 


Vc 


501 ± 44 


2.970 ± 0.005 


15.3 


XcO 


230 ± 40 


3.406 ± 0.007 


6.3 


Vc(2S) 


311 ±42 


3.626 ± 0.005 


8.1 


X(3940) 


266 ± 63 


3.936 ±0.014 


5.0 



The new state has a mass that is above both the DD 
and D*D thresholds. We therefore perform a search 
for A(3940) decays into DD and D*D final states. Be- 
cause of the small product of D^*' reconstruction efficien- 
cies and branching fractions, it is not feasible to recon- 
struct fully the chain e+e" -> J/tp A(3940), A(3940) -> 
D(*W. To increase the efficiency, only one D meson 



in the event is reconstructed in addition to the recon- 
structed J/tp and the other D or D* is detected as a 
peak in the spectrum of masses recoiling against the 
J /4>D combination. The Monte Carlo simulation for 
e + e~ — * J/t/jDD and e + e~ — > J/ipD*D processes in- 
dicates a M ieco i\(J/ijjD) resolution of about 30MeV/c 2 
and a separation between these two processes of 2.5 a. 
Figure [3 shows the M leco ii(J/ipD) spectrum in the D 
mass window (points with error bars) and in the scaled 
D mass sidebands (hatched histogram), where D includes 
D° and D + . Some events have multiple D candidates. 
In these cases, only the candidate with invariant mass 
closest to the nominal D-meson mass is used. Two peaks 
around the nominal D and D* masses are clearly visible 
in this distribution. The excess of real D events com- 
pared to the D sidebands at masses above 2.1 GeV/c 2 is 
due to e+e" -> J/ipD*D* or e+e" -> J/ipD^D^n pro- 
cesses. A fit to this spectrum is performed using shapes 
fixed from MC for three processes (J/ipDD, Jji\)D*D 
and J/ipD*D*) and a second order polynomial. The fit 




FIG. 2: The distribution of masses recoiling against the recon- 
structed J/ipD combinations in the data. Points with error 
bars correspond to the D mass signal window; hatched his- 
tograms show the scaled D sideband contribution. The solid 
line represents the fit described in the text. The dashed line 
shows the background function. 

gives N D p = 86 ±17 (5.1 a) and N D ,^ = 55 ±18 (3.3 a) 
events in the D and the D* peaks, respectively. Select- 
ing events from the M reco n(J/ifjD) regions around the D 
and D* masses (±70MeV/c 2 ), we thus effectively tag the 
processes e + e~ — > J/i/jDD and e + e~ — > J/4>D*D. 

We constrain M reC oii( J/ipD) to the D^ nominal mass, 
thereby improving the M(D^D) = M Teco n(J /ip) resolu- 
tion by a factor of 2.5 (a ~ 10MeV/c 2 after constraint), 
according to the MC simulation. In the A" (3940) -> D*D 
case, the reconstructed D can be from either the A (3940) 
decays or the D* decay: the constraint M roC oii( J/tpD) — > 
M(D*) also works in the latter case, as both A(3940) — > 
D*D and D* decays have very little available phase 
space. 

The resulting M Teco n(J / ip) distributions are shown in 
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Figs. E|a) (M(D) region) and 0b) (M(D*) region). 
The cross-hatched histograms show the scaled D side- 
band distributions. For events with multiple entries, 
the candidate with invariant mass closest to the nomi- 
nal D-meson mass (for D signal window) or closest to 
the center of the sideband (for sidebands) is used. An 
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FIG. 3: The M xeco n(J/tp) distribution for events tagged and 
constrained as a) e + e~ — > J/ipDD, and b) e + e~ — > J/ipD* D. 
The hatched histograms correspond to scaled D sidebands. 
The solid lines are result of the fits, described in the text. 
The dashed lines show: a) the 90% C.L. upper limit on the 
signal; b) the background function. 

A (3940) peak with a resolution that is better than that 
for the unconstrained M YOCO n(J / ip) distribution is evi- 
dent in Fig.|3]b), corresponding to the decay A (3940) — > 
D*D. We perform a fit to this distribution. The sig- 
nal function is a convolution of a Breit-Wigner with a 
free width and a resolution function fixed to the MC ex- 
pectation. The background function is a threshold func- 
tion (A + B ■ M(D*D))^M(D*D) - M thr with M thr = 
M(D*) + M(D). The fit yields the number of signal 
events in the peak N = 24.5 ± 6.9 with a statisti- 
cal significance of 5.0 a. The width of the A(3940) is 
r = (15.4 ± 10.1) MeV/c 2 . The mass of the state is mea- 
sured to be M = (3.943 ± 0.006) GeV/c 2 . 

We perform a similar fit to the M rcco u(J / ip) distri- 
bution for events tagged and constrained as e + e _ — > 
J/ipDD. Since no A(3940) signal is seen for this mode, 
we fit this distribution with A (3940) parameters fixed to 
the values found by the fit of tagged e + e~ — > J/ipD*D. 
The signal yield is found to be O.2+35 events and we set 
an upper limit for the X(3940) signal of 8.1 events at the 
90% C.L. 

An enhancement with a similar mass, F(3940), decay- 
ing into J/ip u> has been recently observed by Belle 0] in 
B decays. We perform a search for the decay A (3940) — > 
J/ipto to see if A (3940) and y(3940) could be the same 



particle. To increase the efficiency we reconstruct the oj 
and only one J/ip from the J/ipJ/tpu final state. The 
unreconstructed J/ip is identified as a peak in the spec- 
trum of recoil masses against the reconstructed J/ipui 
combinations. A signal for A(3940) — > J/ipoj would be 
seen as a peak around the nominal A (3940) mass in a 
distribution of M IBCO a(J/ip) - M Iecoi i(J/ip(j) + M(J/ip) 
if the reconstructed J/ip is prompt, and in the J/ipui 
invariant mass distribution if the reconstructed J/ip is 
from the A (3940) decay. Since the first case has much 
larger combinatorial background and less sensitivity, we 
use only the second case. A scatterplot of M ieco u(J / ipu>) 
vs. J/ip ui invariant mass in the data is shown in Fig.0]a). 
An M{J/ipu>) projection with the additional require- 
ment \M iecoi i(J/ipLu) - Mj/^\ < 100MeV/c 2 is shown 
in Fig. 0]b). A fit to this distribution is done with the 
signal function and parameters fixed from the result of 
the D*D tagged fit; the background is a threshold func- 
tion. The fit yields 1.9I.2' 2 signal events corresponding 
to a 7.4 event upper limit at the 90% C.L. 
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FIG. 4: a) The scatterplot of M rocoi i( J/tp u) vs. M(J/ipu), 
and projection onto b) M(J/ipu>). Points with errors bars 
show the contribution horn the u> mass window; the hatched 
histogram shows the ui sideband. The solid line represents 
the fit described in the text and the dashed line the 90% C.L. 
upper limit on the X(3940) — ► J/ipu contribution. 

The systematic errors for the e + e~ — > J/V'A(3940) 
Born cross section and for the A (3940) branching frac- 
tions £>(A(3940)) are summarized in Table |nj To es- 
timate the systematic errors associated with the fitting 
procedure we study the difference in A (3940) yield re- 
turned by the fit to the M reco n(J/ip) distribution un- 
der different assumptions for the signal and background 
parameterization. In particular, in the first fit we use 
a background function that includes several threshold 
functions corresponding to the production of D*D and 
D*D*. We also use the threshold function (A + B ■ 
-^recoil (J/V'))\/M-ccoii(J/V ; ) - Mthr- Different angular 
distributions result in different J/ip (and D) reconstruc- 
tion efficiencies. In the MC the J/ip production angle and 
J/ip, A (3940) helicity angle distributions are assumed 
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TABLE II: Contribution to the systematic error for 
cr B om(e + e- -> J/^X(3940)) and S(X(3940)) [%]. 



Source 


OBorn 


2?(X(3940)) 


D*D 


DD 


J/lpUJ 


Fitting procedure 


±11 


±17 








Angular distributions 


±19 


±12 




±12 


±16 


A^ch requirement 


±3 


±3 




±3 


±3 


Reconstruction 


±2 


±6 




±6 


±5 


Identification 


±3 


±1 




±1 




Total 


±23 


±22 


±14 


±17 



to be flat. The possible extreme angular distributions 
(1 ± cos 2 8 and sin 2 9) are considered to estimate the sys- 
tematic uncertainty of this assumption. This uncertainty 
partially cancels out in the determination of £?(X(3940)) 
because of the common J/if) efficiency. Other contribu- 
tions come from > 4 requirement efficiency, track 
reconstruction efficiency, lepton identification for recon- 
structed J/ip and kaon identification for reconstructed 
D. 

The systematic errors in the measurement of the 
AT (3940) mass are dominated by the 5 MeV/c 2 uncer- 
tainty associated with the fitting procedure. The un- 
certainty due to the J/ip momentum scale is less than 
3 MeV/c 2 0. These contributions are added in quadra- 
ture to give 6 MeV/c 2 total. The systematic error on the 
mass of the nearby r] c (2S), which is still poorly known ||, 
is found to be the same. From the fit to Fig. |3]b) we es- 
timate the AT (3940) width to be smaller than 47 MeV/c 2 
at the 90% C.L.; this takes into account the fact that the 
likelihood function is not parabolic. When fitting sys- 
tematics are taken into account, we find T < 52 MeV/c 2 
at the 90% C.L. 

The Born cross section for e + e~ — > J/ifjX(39A0) is 
calculated following the procedure used in Ref. |7j. As 
in the previous Belle papers, because of selection criteria 
the result is presented in terms of the product of the 
cross section and the branching fraction of the AT (3940) 
into more than two charged tracks (<8>2)- We obtain 

cTBom x B >2 = (10.6 ±2.5 ±2.4) fb. (2) 

Using the X(3940) yields in inclusive and D*D tagged 
samples, we calculate the fraction of AT(3940) decays with 
more than two charged tracks in the final state into D*D, 
6 >2 (A'(3940) -> D*D). To remove the correlation be- 
tween the two samples, we apply a veto on D* D tagging 
in the inclusive sample. Correcting for the tagging and 
veto efficiencies obtained from MC with equal fractions 
of AT(3940) -> D*°D° and AT(3940) -> D*+D~ assumed, 
we find 

£ >2 (A:(3940) -» D*D) = (96±i ± 22)% 

(> 45% at 90% C.L.), (3) 



where the systematic errors are taken into account for 
the lower limit. In the limit of a vanishing fraction of 
low charged multiplicity X(3940) decays, the measured 
value of B >2 corresponds to S(X(3940) -► D*D). 

We set upper limits on the branching fractions of decay 
of X(3940) into DD and X(3940) -> J/tpuj final states, 
taking into account the estimated systematic errors: 

£(A:(3940) -» DD) < 41% at 90% C.L.; (4) 
£(AT(3940) -> J/tpuj) < 26% at 90% C.L. (5) 

These limits assume that low charged multiplicity 
AT (3940) decays are negligible and, thus, may be over- 
estimated. 

In summary, we have observed a new charmonium- 
like state, AT (3940), produced in the process e + e~ — > 
J I V> X (3940), both in inclusive production and via the 
AT (3940) — ► D*D decay mode. Both observations have a 
5a statistical significance. We have measured the Born 
cross section for the production process, the branching 
fraction for AT(3940) — > D*D, and set upper limits on 
X(3940) decays to DD and J/tpu. 
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